Constraint on Coupled Dark Energy Models from Observations 



o 
o 

> 

o 
in 

(N 

O 

u 

43 : 

Oh- 
i 

o 

o3 



> 
O 
(N 
00 

^1- 



o 



13 



J un- Qing XisQ 

Scuola Internazionale Superiore di Studi Avanzati, Via Beirut 2-4, 1-340 14 Trieste, 

The coupled dark energy models, in which the quintessence scalar field nontrivially couples to 
the cold dark matter, have been proposed to explain the coincidence problem. In this paper we 
study the perturbations of coupled dark energy models and the effects of this interaction on the 
current observations. Here, we pay particular attention to its imprint on the late-time Integrated 
Sachs-Wolfe (ISW) effect. We perform a global analysis of the constraints on this interaction from 
the current observational data. Considering the typical exponential form as the interaction form, 
we obtain that the strength of interaction between dark sectors is constrained as /3 < 0.085 at 
95% confidence level. Furthermore, we find that future measurements with smaller error bars could 
improve the constraint on the strength of coupling by a factor two, when compared to the present 
constraints. 

PACS numbers: 98.80.Cq, 98.80.-k 



I. INTRODUCTION 

Current cosmological observations, such as the cosmic 
microwave background (CMB) measurements of temper- 
ature anisotropies and polarization [l[ and the redshift- 
distance measurements of Type la Supernovae (SNIa) at 
z < 2 I, have demonstrated that the Universe is now 
undergoing an accelerated phase of expansion and that 
its total energy budget is dominated by the dark energy 
component. The nature of dark energy is one of the 
biggest unsolved problems in modern physics and has 
been extensively investigated in recent years, both under 
the theoretical and the observational point of view. 

The simplest candidate of dark energy is the cosmolog- 
ical constant, where the equation of state w always re- 
mains -1. In the standard A-Cold Dark Matter (ACDM) 
cosmology, cold dark matter only interacts with other 
components by gravity, while dark energy is simply a 
cosmological constant without any evolution and pertur- 
bation. Although this concordance model can fit the cur- 
rent observational data very well [J, Q , the possibility of 
dynamical dark energy models cannot be ruled out yet. 
Furthermore, this model suffers of the fine-tuning and co- 
incidence problems, i.e. why the Universe is dominated 
by dark energy in late times [!, H[ . 

In order to lift these severe tensions, many alterna- 
tive dynamical dark energy models, such as q uintessence 
@, @JZL 1, phantom 0], K -essence [lfj, [ll| and quin- 
tom 1121 . have been proposed recently (for a review see 
Ref. [131] V For example, the quintessence model, a scalar 
field <j> slowly rolling down a potential energy V ((/>), has 
the spatial fluctuations and its equation of state can 
evolve with the cosmic time. More importantly, being a 
dynamical component, the quintessence is naturally ex- 
pected to interact with the other components, such as the 
cold dark matter [3, [l5[ or massive neutrinos [H, [I?) , 
in the field theory framework. If these interactions really 



exist, it would open up the possibility of detecting the 
dark energy non-gravitationally. 

In the coupled dark energy models, the quintessence 
scalar field could nontrivially couple to the cold dark mat- 
ter component. The presence of the interaction clearly 
modifies the cosmological background evolutions. The 
evolution of cold dark matter energy density p c (a) will 
be different, when compared to that of the minimally 
coupled models, and dependent on the quintessence field 



p e (o) = Pc0 a- S+S ^ 



(1) 
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where a is the scale factor, p cQ is the present value of cold 
dark matter energy density, and 5(<p) is the modification 
due to the interactions and could be non-zero during the 
evolution of Universe. In this case, at early times there 
will be more (less) cold dark matter energy density, when 
S{<f>) <0{>0). 

On the other hand, the interaction between dark sec- 
tors will also affect the evolution of cosmological pertur- 
bations, which has been widely investigated recently (see 

e.g. Refs.Ql m M M H M MM M M El). 

Due to the different evolution of cold dark matter en- 
ergy density, the interaction between dark sectors could 
shift the matter-radiation equality scale factor a eq , and 
affect the locations and amplitudes of acoustic peaks of 
CMB temperature anisotropies and the turnover scales of 
large scale structure (LSS) matter power spectrum con- 
sequently. In addition, the interaction will also affect 
the late ISW effect [29| at large scales which is produced 
by the CMB photons passing through the time-evolving 
gravitational potential well, when dark energy or curva- 
ture becomes important at later times (23l . [30| . 

Therefore, with the accumulation of observational data 
and the improvements of the data quality, it is of great 
interest to investigate the non-minimally coupled dark 
energy models from the current observational data. In 
this paper we study the perturbations of coupled dark 
energy models and present the constraints on the inter- 
actions from the observational data in detail. The struc- 
ture of the paper is as follows: in Sec. II and Sec. Ill we 
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show the basic equations of background evolution and 
linear perturbations of the coupled dark energy model, 
respectively. In Sec. IV we present the current and fu- 
ture observational datasets we used. Sec.V contains our 
main global fitting results from the current and future 
observations, while Sec. VI is dedicated to the summary. 



II. BACKGROUND EVOLUTION 

In our analysis we assume a flat, homogeneous, 
Fricdmann-Robcrtson- Walker universe with metric 



ds 2 = a 2 (77) {—drj 2 4- dxidx 1 



(2) 



where rj is the conformal time. The Friedmann equation 
is: 



8ttG 



(3) 



where p tot is the total energy density of Universe and the 
dot refers to the derivative with respect to the conformal 
time 77. 

From the Lagrangian of the quintessence scalar field: 



we can obtain the energy density and pressure 

And the Friedmann equation Eq.((3]) becomes: 
BttG , / 



H 2 = 



Pb 



(4) 



(5) 



(6) 



where p 7 , p h and p c are the energy densities of radia- 
tion, baryons and cold dark matter, respectively. Con- 
sequently, we define the energy density parameters f2 ; = 
Pi/ptot, where p\ is the energy density of each compo- 
nent. In the following calculations we will set the reduced 
Planck mass M p i = 1/ V&nG = 1 . 

Here, we consider the typical exponential form as the 
quintessence potential, namely 



V{4>) = V e 



-\<t> 



We also take the exponential form: 



(7) 



(8) 



as the interaction form between quintessence and cold 
dark matter, where p* is the bare energy density of cold 
dark matter and (3 is the strength of interaction. There 
are several stringent constrains on this strength of inter- 
action, namely, \(3\ < 0(0.1) at 95% confidence level, 
from different observations 0, [3l|, H^. Using these 



potential forms, the initial choice of <f>i will not affect 
the evolution of Universe significantly. And this cou- 
pled model could be reduced to the standard ACDM 
model by choosing A = (3 — straightforwardly. The 
effective potential of this coupled system is given by: 
V e g(4>) — V((f>) + p c {4>)- If we restrict to the case 
A > 0, this effective potential will have a minimal value at 
0min = hi (XVo//3p c )/X for the strength (3 > 0. Therefore, 
in our following calculations we will restrict the forms to 
the case A > and (3 > to lead to the acceleration 
expansion of Universe at late times. 

When including the interactions, the conservation of 
energy momentum for each component becomes p3l |: 



(9) 



In our analysis, we assume that the baryons and radiation 
are not coupled with the quintessence scalar field. There- 
fore, when considering the interaction between cold dark 
matter and quintessence, the energy conservation equa- 
tions of CDM and quintessence will be violated, while 
those of baryons and radiation are held: 

p 7 + m Pl = , (10) 

Ph + 3Hp h - , (11) 
Pc + m Pc = Q c = {3j> Pc , (12) 
p + 3Hp^l + w+) = Q<t, = -[3<j>p c , (13) 

where the energy exchange Q c is the function of 4>. When 
Q c < (> 0), the energy of cold dark matter (dark en- 
ergy) transfers to dark energy (cold dark matter). Be- 
cause we describe the dark energy component using a 
quintessence scalar field, the energy conservation equa- 
tion Eq. (fT3)) will be describe as the Klein-Gordon equa- 
tion: 



4> + 2H4> + a 2 V'{(f) = -a 2 Pp c 



(14) 



where the prime denotes the derivative with respect to 
the quintessence scalar field <j>: V'{(j>) = dV (<j>) / d(j> . 

We can also define the effective equation of state of 
CDM and quintessence field to describe the equivalent 
uncoupled model in the background: 



Pc 

P<t> - 



-m Pc {\ 
mpAi 



wf) = 0, 

„..eff\ 



') = 0, 

where the effective equation of state are given by: 



w, 



cff 



cff 



W<f, 



P<j> Pc 

3H Pa, 



(15) 
(16) 



(17) 



We can find that in this case the effective equation of 
state of cold dark matter w^r will be non-vanishing and 
evolves with the cosmic time. 

Based on these equations, we can study the back- 
ground evolutions of each component in the non-minimal 
coupling system. In Fig[T]we illustrate the evolutions of 
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FIG. 1: The cosmological evolutions of the cold dark mat- 
ter and quintessence energy density parameters, fii, for two 
models: A = 1.22, p = (black solid lines) and A = 1.22, 
(3 = 0.15 (red dashed lines). 



the cold dark matter and quintessence energy density pa- 
rameters as a function of redshift for two different models: 
the uncoupled model A = 1.22, (3 — (black solid lines) 
and the coupled system A = 1.22, (3 — 0.15 (red dashed 
lines). For other cosmological parameters, we fix them to 
the best fit values of WMAP5 data [l|: O b /i 2 = 0.02267, 
il c h 2 — 0.1131, and h = 0.705. Here, we normalize en- 
ergy density parameters to their present values. 

Clearly, the interaction between the quintessence 
scalar field and the cold dark matter significantly modi- 
fies the evolutions of their energy density. Because of the 
presence of the coupling, during the evolution, j3<j> < 
and the energy density of cold dark matter evolves faster 
than a -3 . Therefore, if we normalize energy density pa- 
rameters to their present values, the cold dark matter 
energy density will be larger than that of the uncoupled 
model. The energy of cold dark matter transfers to dark 
energy. 



III. LINEAR PERTURBATION 

In this section we will consider the evolution of cos- 
mological perturbations in the non-minimally coupled 
quintessence model. Starting from the metric in the syn- 
chronous gauge [H, [H, [3f| , the line element is given by: 



ds 2 



-a 2 (rj) (-drj 2 + (<% + hij)dx l dx j ) 



(18) 



where hij is the metric perturbation in this synchronous 
coordinate system. Here, we restrict only for the scalar 
mode of the metric perturbations. 

We write the inhomogcncous energy density of cold 
dark matter and the quintessence scalar field as: 

P c (v,x) = p c (rj) + Sp c (r),x) 



(j)(r),x) 



(19) 
(20) 



where p c (r]) and <j){rj) are the background parts, while Sp c 
and S(f> are the perturbations. Using the perturbed part 
of the energy momentum conservation equation Eq.((5]) 
for the coupled system, we could calculate the evolution 
equations for the perturbation of cold dark matter: 



Spc 



S c = 3{H +[3(t>)[w c --^)6 c -(l 



+f3{l - 3w c )6(/> + 04>5(j>{\ - 3w c ) , 

w c a . Sp c /5p c 




= -«(l-3u) c )0 c 



1 



1 



k 2 S c 



w c 



1 + W c 



k 2 S<j) - (3(1 - 3w c 



(21) 



(22) 



where Sp c is the pressure perturbation of cold dark mat- 
ter, h is the metric perturbation, 9 C = ik^v^ is the gra- 
dient of the velocity field, and a c is the shear stress of 
cold dark matter. In our coupling form Eq.([5]), (3 is a 
constant, so we have (3' = d/3/d<f> = 0. Finally we obtain 
the perturbation equations of cold dark matter: 



5 C = -9 c -h/2 



: + k 2 (35(() 



(23) 
(24) 



We note that in the presence of interaction, the gradi- 
ent of the velocity 9 C will evolve to be nonzero, even if 
its initial value is zero; whereas 9 C remains zero at all 
times in the minimally coupled model. Therefore, in our 
calculations we will evolve the perturbation equations in 
an arbitrary synchronous gauge, instead of the cold dark 
matter rest frame [24J . 

Meanwhile, the perturbed Klein-Gordon equation is 
given by: 

64> + 2?H5j) + k 2 S(j) + a 2 V"S(j) + h(f,/2 = -a 2 (3p c 8 c . (25) 

Wc find that the perturbations of this non-minimally cou- 
pled system are stable. This is different from some cou- 
pling models in which dark energy component is modelled 
as a fluid with constant equation of state parameter w. 
In those models, when w is close to the cosmological con- 
stant boundary, the coupling term Q ~ p c will lead to an 
instability [H, E3, El Hi| . 

In this paper we still choose the adiabatic initial con- 
dition, which implies that the entropy perturbation be- 
tween cold dark matter and quintessence: 



5=^-^=0 
Pc Pa, 



(26) 



and the quintessence intrinsic entropy perturbation: 



(27) 
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FIG. 2: The CMB temperature anisotropies and LSS matter power spectrum for two models: A = 1.22, = (black solid 
lines) and A = 1.22, = 0.15 (red dashed lines). 



vanish at early times. However, during the evolution 
of Universe, the iso-curvature perturbations will be pro- 
duced and the adiabatic condition is generally no longer 
conserved, even for large-scale modes, due to the presence 
of non- minimal coupling [2ol |28| . 

Following the above equations, we modify the CAMB 
code [33] to calculate the CMB temperature power spec- 
trum Cj T and LSS matter power spectrum P(k) in 
this non-minimal coupling scenario. Besides the back- 
ground parameters, we also set r = 0.084, n s = 0.96 and 
A s = 2.15 x 10~ 9 at k = 0.05 Mpc -1 , which are the best 
fit values obtained from the WMAP5 data [l|. 

In Fig[2]we plot the CMB temperature power spectrum 
and LSS matter power spectrum for two different models. 
In this non-minimally coupled system the matter power 
spectrum is almost unchanged and the growth of density 
perturbations are not affected significantly on the large 
scales [2l|. However, on the small scales the situation 
will be different. As we mentioned before, at early times 
there is higher cold dark matter energy density in the 
coupled model with j3 > 0. In this case, the epoch of 
equality a oq occurs further from the recombination. An 
earlier a oq means only the very small scale modes enter 
the horizon and decay during the radiation dominated 
epoch. Therefore, the amplitude of small scale matter 
power spectrum is enhanced and the turnover occurs on 
the smaller scales as shown in the left panel of FigfS] 
Consequently, the expected value of a% will become larger 
than that of the minimally coupled case. 

In the CMB temperature power spectrum, the most 
obvious effect is that the amplitude on the small scales 
has been decreased significantly due to the presence of 
interaction between dark sectors. In the non-minimally 
coupled model, the epoch of equality a cq occurs further 
from the recombination, so that at recombination the in- 
homogeneities induced from the radiation becomes small. 
The small-scale anisotropies decrease when the coupling 



strength is increasing. In addition, the locations of the 
acoustic peaks in the CMB temperature anisotropies will 
also slightly shifted to the smaller scales, due to the pres- 
ence of interaction. 

More interestingly, the non-minimal coupling [3 > 
also suppresses the anisotropies on the large scales {I < 
20) which arise mainly from the late-time ISW effect. 
When /3 is positive, the energy transfers from cold dark 
matter to dark energy and at early times cold dark matter 
density becomes larger. Consequently, the dark matter- 
dark energy equality scale factor a' eq will be larger than 
that of the un-coupled case. Because the gravitational 
potential $ stays constant in the matter dominated era, 
$ and additional CMB anisotropies at large scales will be 
suppressed in the presence of non-minimal coupling. The 
CMB temperature anisotropies on very large scales be- 
come small, as shown in Figl2] Therefore, observing the 
late time ISW effect could be a helpful way of studying 
the nontrivially coupling between dark sectors. 

However, the most significant ISW effect contributes 
to the CMB anisotropies on large scales that are strongly 
affected by the cosmic variance. Fortunately, this prob- 
lem can be solved by the cross-correlation between ISW 
temperature fluctuation and the density of astrophysi- 
cal objects like galaxies or quasars. (For the details, we 
refer the reader to Refs.jH, [H, |4(j.) In Fig [3] we plot 
the cross-correlation power spectra Cf T for two different 
models using the public package CAMB-sources 1 . One 
can find that the power spectrum Cf T of the non-minimal 
coupling case will also be suppressed. The non-minimal 
coupling leaves an interesting imprint on the cross corre- 
lation power spectrum. We expect that the observational 
ISW data could improve the constraints on the interac- 
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only determine an effective distance [49(: 
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FIG. 3: The cross-correlation signal C e for two models: A = 
1.22, (3 = (black solid lines) and A = 1.22, (3 = 0.15 (red 
dashed lines). 



tion. 



IV. OBSERVATIONAL DATA 



A. Current Datasets 

In our calculations, we will rely here on the following 
cosmological probes: i) CMB anisotropies and polariza- 
tion; ii) baryonic acoustic oscillations in the galaxy power 
spectra; iii) SNIa distance moduli; and iv) the angu- 
lar auto-correlation (ACF) and cross-correlation (CCF) 
functions data obtained from the quasar map and CMB 
map. 

In the computation of CMB power spectra we have 
included the WMAP five-year (WMAP5) temperature 
and polarization power spectra with the routines for 
computing the likelihood supplied by the WMAP team 
[2, EH El, EE 0, IH. Since the non- minimal cou- 
pling also leaves some signatures on the small scales 
CMB power spectrum, in our calculations we also in- 
clude some small-scale CMB measurements to improve 
the constraints, such as BOOMERanG 0, CBI [43] and 
ACBAR Jig. 

BAOs (Baryonic Acoustic Oscillations) have been de- 
tected in the current galaxy redshift survey data from the 
SDSS and the Two-de gree Field Galaxy Redshift Survey 
(2dFGRS) [H, M, HI 113 • The BA0 can directly mea- 
sure not only the angular diameter distance, Da(z), but 
also the expansion rate of the universe, H(z), which is 
powerful for studying dark energy [13 ■ Since current 
BAO data are not accurate enough for extracting the in- 
formation of Da (z) and H{z) separately [13, one can 



(l + zfD\(z) 



H(z) 



1/3 



(28) 



In this paper we use the Gaussian priors on the distance 
ratios r s {z d )/D v (z): 



r s (z d )/D v (z 
r s (z d )/D v (z 



0.20) = 0.1980 ±0.0058 
0.35) = 0.1094 ±0.0033 



(29) 



with a correlation coefficient of 0.39, extracted from the 
SDSS and 2dFGRS surveys [13, where r s is the comoving 
sound horizon size and z d is drag epoch at which baryons 
were released from photons given by Ref. [H[ . 

The SNIa data provide the luminosity distance as a 
function of redshift which is also a very powerful measure- 
ment of cosmological evolution. The supernovae data we 
use in this paper are the recently released Union com- 
pilation (307 samples) from the Supernova Cosmology 
project [2J, which include the recent samples of SNIa from 
the (Supernovae Legacy Survey) SNLS and ESSENCE 
survey, as well as some older data sets, and span the 
redshift range < z < 1.55. In the calculation of the 
likelihood from SNIa we have marginalized over the nui- 
sance parameter as done in Ref. 56]. 

For the observational ISW data, we use the ACF and 
CCF data obtained from the quasar map of Sloan Dig- 
ital Sky Survey Data Release Six H3 and the WMAP5 
Internal Linear Combination (LLC) map |4ll|. which has 
been used in Ref. [13- (Also, Refs. [59L l60j] provided the 
observational ISW data using different CMB and LSS 
surveys.) This ISW data set gives the information at 
very low redshift {z ~ 1) which is much lower than the 
recombination [z ~ 1100). Thus, we do not consider the 
possible covariance matrix between this ISW data and 
the WMAP power spectra Cg. 

Furthermore, we make use of the Hubble Space Tele- 
scope (HST) measurement of the Hubble parameter 
Hq = 100 h km s _1 Mpc -1 by a Gaussian likelihood 
function centered around h = 0.72 and with a standard 
deviation a = 0.08 [H. 



B. Future Datasets 

In order to forecast future measurements we will use 
the same observables as before without BAO. 

For the simulation with Planck [13, we follow the 
method given in Ref. [13 and mock the CMB temper- 
ature (TT) and polarization (EE) power spectra and 
temperature-polarization cross-correlation (TE) with the 
isotropic noise by assuming a given fiducial cosmologi- 
cal model. In Table I, we list the assumed experimental 
specifications of the future (mock) Planck measurement. 
Here we neglect foregrounds and the tensor information. 
The effective x 2 is: 



Xcff 



-21og£ 
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TABLE I. Assumed experimental specifications for the mock 
Planck-like measurements. The noise parameters At and 
Ap are given in units of /^K-arcmin. 



/sky 


f 

*rnax 


(GHz) 


^fwhm 


A T 


Ap 


0.65 


2500 


100 


9.5' 


6.8 


10.9 






143 


7.1' 


6.0 


11.4 






217 


5.0' 


13.1 


26.7 



los [ ww^W^f ) (30) 

r*TTr<EE , r'TTriEE or^TEr'TE ~] 
C TT C EE — (C TE ) 2 i ' 

where Cf Y denote theoretical power spectra and C* Y 
denote the power spectra from the simulated data. The 
gaussian likelihood function has been normalized with 
resp ect to the maximum likelihood, where Cf Y = Cf Y 

The proposed satellite SNAP (Supernova / Accelera- 
tion Probe) will be a space based telescope with a one 
square degree field of view that will survey the whole sky 
[66j |. It aims at increasing the discovery rate of SNIa to 
about 2000 per year in the redshift range 0.2 < z < 1.7. 
In this paper we simulate about 2000 SNIa according to 
the forecast distribution of the SNAP j22|. For the er- 
ror, we follow the Ref. [63| which takes the magnitude 
dispersion to be 0.15 and the systematic error cr sys = 
0.02 x z/1.7. The whole error for each data is given by 

CTmagC^i) = J~o%J&) + 0.15 2 /ni , where n 4 is the num- 
ber of supernovae of the i'th redshift bin. Here, we do 
not consider the possible covariance matrix among these 
redshift bins. 

For the future ISW ACF and CCF data, we simulate 
the mock dataset from the best fit values of current data 
combination CMB+BAO+SNIa+ISW. We reduce the er- 
ror bars of these data points by a factor three. We also 
use directly the covariance matrix taken from present 
data and divide it by a factor nine. This improvement 
could be achievable by next generation of large-scale sur- 
veys such as the Large Synoptic Survey Telescope (LSST, 
|68|), the Panoramic Survey Telescope and Rapid Re- 
sponse System (Pan-STARRS, [69() and the Dark Energy 
Survey (DES; The Dark Energy Survey Collaboration 
2005). These surveys are likely to allow for an order of 
magnitude improvement in the number of quasars, which 
account for the factor three considered here if we assume 
that the error bars scale as yj -/V quasar . Also, the coverage 
of the sky fraction is at present of the order of 20% by 
SDSS DR6, so an all-sky survey will already give a factor 
two improvement. Another possible improvement in this 
direction is the use of type-2 quasars instead of the type- 
1 used here that could further decrease the error bars of 
the sample (see the discussion in Ref . (57| ) . 
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FIG. 4: One dimensional posterior distributions of the 
strength of interaction (3 from various data combinations: 
WMAP5 alone (black solid line), WMAP5+CMBExt (red 
dashed line) , WMAP5+ISW (blue dash-dot line) , all current 
data sets (magenta short dashed line), and all future data sets 
(navy blue short dash-dot line). 



V. NUMERICAL RESULTS 

In our analysis, we perform a global fitting using 
the CosmoMC package [70| a Monte Carlo Markov Chain 
(MCMC) code, which has been modified to calculate the 
theoretical ACF and CCF. We assume purely adiabatic 
initial conditions and a flat universe, with no tensor con- 
tribution. Besides the parameters A and (3 in the poten- 
tial and coupling forms, we also vary the following cos- 
mological parameters with top-hat priors: the dark mat- 
ter energy density f2 c /i 2 e [0.01,0.99], the baryon energy 
density VL^h 2 g [0.005, 0.1], the primordial spectral index 
n s € [0.5,1.5], the primordial amplitude log[10 10 A s ] S 
[2.7, 4.0] and the angular diameter of the sound horizon 
at last scattering € [0.5, 10]. For the pivot scale we set 
fc s0 = 0.05 Mpc -1 . When CMB data are included, we 
also vary the optical depth to reionization r G [0.01, 0.8]. 
We do not consider any massive neutrino contribution. 
From the parameters above the MCMC code derives the 
reduced Hubble parameter Hq, the present matter frac- 
tion f2 m o, and as, so, these parameters have non-flat pri- 
ors and the corresponding bounds must be interpreted 
with some care. There are three more parameters re- 
lated to the ACF and CCF data: the constant bias b, 
the efficiency of quasar catalog a and the ISW amplitude 
A amp which is defined by: C qT ((9) = A amp C qT (60, where 
C qT and C qT are the observed and theoretical CCF. In 
addition, CosmoMC imposes a weak prior on the Hubble 
parameter: h £ [0.4,1.0]. 

In Fig|4]we show the one-dimensional posterior distri- 
butions of [3 from various data combinations. Firstly, 
we present the constraint from the WMAP5 data alone 
(black solid line). For the strength of interaction, the 
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constraint is still weak, namely the 95% upper limit is 
(3 < 0.165. which is consistent with other results in the 
literature [24], [3l|, [32j • As we mentioned before, the non- 
minimal coupling shifts the acoustic peaks of CMB tem- 
perature anisotropies on the small scales. The small-scale 
CMB measurements (CMBExt) should improve the con- 
straints on the coupling strength. Therefore, we plot the 
constraint on /3 from WMAP5+CMBExt data combina- 
tion in FigfJ] (red dashed line). When adding the small- 
scale CMB data, we can find that the constraint on j3 
improves slightly: (3 < 0.135 at 95% confidence level. 
The small-scale CMB observations indeed improve the 
constraint. 

We also include the observational ISW data 58] into 
our calculations. In Fig|4]we show the one dimensional 
distribution of (3 from WMAP5+ISW data combination 
(blue dash-dot line). We can find that the combined 
constraint from WMAP5+ISW is slightly improved over 
using WMAP5 alone, namely, the 95% upper limit of 
the coupling strength is (3 < 0.130. Since at present 
constraints from the ISW data are still very weak and 
in the calculations we only consider the high-redshift 
quasar catalog and totally neglect other low-redshift trac- 
ers which could give powerful ISW constraints [59|, H3| , 
the constraint on (3 does not improve significantly. When 
combining WMAP5, CMBExt and ISW data together, 
the constraint on (3 will improve further: (3 < 0.120 at 
95% confidence level. 

Finally we combine all the current datasets together 
to constrain the coupling strength (3. In FigfJ] we find 
that the constraint tightens significantly (magenta short 
dashed line): 

< 0.050 (0.085) (31) 

at 68% (95%) confidence level. The 95% upper limit is re- 
duced by a factor of 2, when compared to the constraint 
from WMAP5 alone, which is due to the constraining 
power of SNIa and BAO on the evolutions of cosmolog- 
ical background parameters. Meanwhile, the parameter 
A in the potential form of quintessence scalar field is con- 
strained to be A < 1.05 at 95% confidence level, which is 
consistent with the result in Ref. [24j . 

As we mentioned before, the non-minimal coupling 
(f3 > 0) affects the large scale structure formation. From 
FigJ5] we can see that the matter power spectrum on the 
small scales will be enhanced significantly by the large 
value of (3. The positive (3 will lead to a high value of a s 
today. Using the all datasets combination, we obtain the 
limit on a s today of a s = 0.816±0.040 (68% C.L.), which 
is obviously higher than one obtained in the minimally 
coupled model (/? = 0): a 8 = 0.793 ± 0.030 (68% C.L.), 
while the error bar is enlarged slightly. Therefore, the as 
today and (3 are correlated as we expected, as shown the 
black solid lines in Fig [5] 

In Fig|6] we also show the two dimensional contour 
in the (A, (3) panel. When A is larger than zero, the 
equation of state of dark energy will evolve with cos- 
mic time, and the dark energy acts more like the dark 



0.95 - 




0.02 0.04 0.06 0.08 0.1 0.12 

P 

FIG. 5: Two dimensional marginalized contour in the (/3, as) 
panel from all current data sets (black solid lines) and all 
future data sets (red dash-dot lines), respectively. 



0.12- 




1.4 



FIG. 6: Two dimensional marginalized contour in the (A, f3) 
panel from all current data sets combination. 



matter (w^ > —1). In this case the dark energy will 
dominate our Universe at relatively higher redshift and 
slow down the linear growth function of the matter per- 
turbation and leads to a low value of <7g today. The con- 
straint of a s = 0.793 ±0.030 (68% C.L.) from all current 
datasets is obviously lower than that of the pure ACDM: 
a s = 0.812±0.026 (68% C.L.) 0. As the parameter A in- 
creases, in order to reach the same value of erg today, the 
strength of coupling (3 should be increased at the same 
time. Therefore, A and (3 are also correlated, illustrated 
in FiglHl 

From the results presented above, we can see the 
current observations could give the constraints on the 
strength of interaction as (3 < 0(0. 1) at 95% confidence 
level. It is worthwhile to discuss whether future data 
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could give more stringent constraints on (3. 

Therefore, in Fig0] we also present the constraint on 
the strength of interaction (3 from the future datasets 
(navy blue short dash-dot line). Due to the smaller error 
bars of the mock data sets, the constraint on (3 improves 
significantly, namely the 68% (95%) upper limits are now 

(3 < 0.025 (0.039) , (32) 

which are reduced by another factor of 2, when com- 
pared to the current constraint Eq. (|3l"j) . We also show 
the two dimensional contour in the er 8 ) panel (red 
dash-dot lines) in FigO The correlation still exists and 
the limits have been shrunk significantly. These results 
imply that the future measurements could constrain the 
non-minimal coupled dark energy model with a higher 
precision. 

VI. SUMMARY 



cold dark matter, could affect the CMB temperature 
anisotropies and LSS matter power spectrum. In this 
paper we present constraints on this coupled dark en- 
ergy model using the latest observations and future mea- 
surements. The WMAP5 data alone could only give a 
weak constraint on the strength of coupling (3. And then 
we exploit the capabilities of the late-time ISW effect in 
constraining the non-minimal coupling, using the cross- 
correlation signal between the quasar sample of SDSS 
DR6 H3 and the WMAP5 ILC map 0. We find that 
the current ISW data could slightly improve the con- 
straint on (3. If we add the BAO and SNIa data into 
our calculations, the constraint on (3 will improves sig- 
nificantly, namely the 95% upper limit is (3 < 0.085. Fi- 
nally we simulate the future measurements with smaller 
error bars and find that the future measurements could 
improve the constraint on the strength of coupling by a 
factor of two, when compared to the present constraints. 



The coupled dark energy models, in which the 
quintessence scalar field non-minimally couples to the 
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